ON to-DIMENSIONAL TORIC CODES 



JOHN LITTLE AND RYAN SCHWARZ 



Abstract. Toric codes are a class of m-dimensional cyclic codes introduced 
recently by J. Hansen in 0, [8], and studied in [§], I1UI . They may be 
defined as evaluation codes obtained from monomials corresponding to integer 
lattice points in an integral convex polytope P C R m . As such, they are in a 
sense a natural extension of Reed-Solomon codes. Several articles cited above 
use intersection theory on toric surfaces to derive bounds on the minimum 
distance of some toric codes with m = 2. In this paper, we will provide a 
more elementary approach that applies equally well to many toric codes for 
all m > 2. Our methods are based on a sort of multivariate generalization 
of Vandermonde determinants that has also been used in the study of multi- 
variate polynomial interpolation. We use these Vandermonde determinants to 
determine the minimum distance of toric codes from rectangular polytopes and 
simplices. We also prove a general result showing that if there is a unimodular 
integer affine transformation taking one polytope Pi to a second polytope P2, 
then the corresponding toric codes are monomially equivalent (hence have the 
same parameters). We use this to begin a classification of two-dimensional 
cyclic toric codes with small dimension. 



1. Introduction 

In 0, J. Hansen introduced the notion of a toric code. Let P C W n be an 
integral convex polytope (the convex hull of some set of integer lattice points). We 
suggest as a good general reference for the geometry of polytopes. Suppose that 
P n Z m is properly contained in the rectangular box [0, q — 2] m (which we denote 
□ g _i), for some prime power q. Then a toric code is obtained by evaluating linear 
combinations of the monomials with exponent vector in P n Z m at some subset 
(usually all) of the points of (F*) m . We formalize this in the following definition. 

Definition 1.1. Let ¥ q be a finite field with primitive element a. For f G 27™ 
with < ft < q - 2 for all i, let pf = (a fl , . . . ,a im ) in (F*) m . For each e = 
(ei, . . . , e m ) G Pfl 17"" , let x e be the corresponding monomial and write 

{p f f = (a* )«*-.. (a' m ) em - 

The toric code Cp(¥ q ) over the field ¥ q associated to P is the linear code of block 
length n = (q — l) m with generator matrix 

G=(( P /) e ), 

where the rows are indexed by the efPfl Z m , and the columns are indexed by the 
Pf G (¥*) m . In other words, letting L — Span{x e : e G P H Z m } ; we define the 
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evaluation mapping 

cv:L -> F^- 1 )™ 

9 - (g(p f ) ■■ f 6 (F ? *) m ) 

TTien Cp = ev(L). If the field is clear from the context, we will often omit it in the 
notation and simply write Cp. The matrix G will be called the standard generator 
matrix for the toric code. 

Because of the close connection between integral polytopes and the theory of toric 
varieties, Hansen and others have proposed techniques from algebraic geometry 
such as intersection theory on algebraic surfaces and higher dimensional varieties 
to study toric codes and their parameters. The articles [7], and EH have 

used this approach. 

In this article, for the most part, we use a more elementary viewpoint, based on 
the observation that the square submatrices of the standard generator matrix of a 
toric code are examples of a multivariate generalization of the familiar univariate 
Vandermonde matrices. These multivariate Vandermonde matrices have been stud- 
ied by a number of different techniques in the context of multivariate polynomial 
interpolation. The literature there is truly vast because of the many applications 
of interpolation in numerical analysis and other parts of applied mathematics. We 
direct the reader to the bibliography of 6 . To our knowledge, these techniques 
have not been used before in coding theory in this form. 

Our contributions are as follows. In §2, we begin with a lower bound on the 
minimum distance of toric codes based on Vandermonde determinants (Proposition 
2.1). We use this to study the minimum distances of toric codes from rectangular 
polytopes and simplices (see Theorems 2.4, 2.5, 2.9). We do this by identifying 
special configurations of points in (F*) m for which the Vandermonde determinant is 
nonzero. In the context of polynomial interpolation, such sets are called poised sets 
for the interpolation problem using linear combinations of monomials corresponding 
to the lattice points in some polytope. (They are the sets for which the interpolation 
problem has a unique solution for all function values assigned at those points.) Our 
methods here are suggested by the interpolation-theoretic computations in 0]. 

All of the 2-dimensional examples from jj] and [B] can be handled with our 
methods. But in fact our theorems are more general since they apply for all m > 2, 
not only the case m — 2. 

In §3 we prove the general statement that lattice equivalent polytopes P\,Pi 
yield monomially equivalent toric codes Cp 1 , Cp 2 . We apply this result to consider 
a classification of m = 2 toric codes up to monomial equivalence for small k. 

Finally, some comments about the utility of toric codes are probably in order. In 
the case m — 1, a toric code is just a Reed-Solomon code since P is a line segment 
in [0, q — 2] C M with integer endpoints. Higher dimensional toric codes are in a 
sense a natural extension of Reed-Solomon codes and have many similar properties. 
For instance, it is easy to see that they are all m-dimensional cyclic codes (see 0). 
So one might hope that toric codes exist having similarly good parameters. And 
indeed, [§] contains a number of examples showing that toric codes can have very 
good parameters, equaling or bettering the best known minimum distance for a 
given n, fc in Pj. Not all toric codes perform this well, however. Our main results 
on toric codes from rectangular polytopes and simplices show, in fact, that their 
minimum distances are often quite small for their dimensions. It is an interesting 
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problem, we believe, to determine criteria concerning which polytopcs yield good 
toric codes. 

We will use this notation in the following sections. Suppose that P C dg-i C K m 
is an integral convex polytope. We will write #(P) for the number of integer lattice 
points in P (that is, #(P) = \P n Z m |). We will write 

PnZ m = {e(i):i = l,...,#(P)} 

for the set of those integer lattice points. For any set ScR", conv(S') denotes the 
convex hull of S. 

2. Minimum distances via Vandermonde matrices 

We begin by describing the Vandermonde matrices involved here. Using the 
notation introduced in §1, let P be an integral convex polytope, and suppose PC] 
7L m = {e(i) : i = 1, . . . , #(P)}, listed in some particular order. Let S = {pj : j = 
1, . . . , #(P)} be any set of #(P) points in (F*) m , also ordered. Then V(P; S), the 
Vandermonde matrix associated to P and S, is the #(P) x #(P) matrix 

V(P;S)=(pf), 

e(i) 

where we use the standard multi- index notation p ■ to indicate the value of the 
monomial x e ^ at the point pj. For example, if P = conv{(0, 0), (2,0), (0,2)} in 
R 2 , and S = {{xj,yj)} is any set of 6 points in (F*) 2 , for one particular choice of 
ordering of the lattice points in P, we have 



(1) V(P;S) = 
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Since we assume P C O q —i, the monomials x e ^ define linearly independent 
functions on (F*) m and we can also view the V(P; S) as square submatrices of the 
standard generator matrix for the toric code Cp = Cp(¥ q ). Our first observation 
is that Vandermonde determinants may be used to bound the minimum distance 
of a toric code. 

Proposition 2.1. Let P C M. m be an integral convex polytope. Let d be a positive 
integer and assume that in every set T C (F*) m with \T\ = (q — l) m — d + 1 there 
exists some S C T with \S\ = #(P) such that det V(P; S) 7^ 0. Then the minimum 
distance satisfies d(Cp) > d. 

Proof. All codewords of Cp are linear combinations of the rows of the #(P) x (q — 
\) m standard generator matrix. If a codeword has zeroes in the locations corre- 
sponding to the subset T c (F*) m , then by looking at the entries corresponding to 
ScTwe get a system of #(P) homogeneous linear equations in the coefficients 
of the linear combination, whose matrix is V(P; S). By hypothesis, this matrix is 
nonsingular, so all the coefficients in the linear combination must be zero. Since 
this is true for all T, all nonzero codewords of Cp have at most (q — l) m — d zero 
entries, which implies d(Cp) > d. □ 
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We will apply this proposition by identifying particular configurations of points 
S for which det V(P; S) ^ 0, based on the particular monomials appearing in 
P. From examples such as (1) above, it should be relatively clear that giving 
characterizations of S such that det V(P;S) — or det V(P;S) ^ is difficult in 
general. Indeed, even in the context of multivariate polynomial interpolation, only 
partial results in special cases are well understood. 

We will begin with the case where P is a rectangular polytope in R m (special 
cases are rectangles in the plane and rectangular solids in R 3 ). For these polytopes, 
the extension from the case m — 2 to general m > 2 is almost immediate, so for 
notational simplicity, we will treat only the case m = 2 in detail. 

Let P k ,e be the rectangle P k . e = conv{(0, 0), (fc, 0), (0, £), (fc, £)}. Note that 
#(Pk,e) = {k + l){£ + I). We will call any set S of (k + 1)(£ + 1) points in (F*) 2 
consisting of (£ + 1) distinct points on each of (fc + 1) distinct vertical lines x = en 
a (fc + 1) x (1 + 1) configuration. 

When we construct V{Pk,t] S) f° r a (fc + 1) x (^+1) configuration we get a matrix 
as in the following proposition. 



Proposition 2.2. Suppose A = (cy) is an a x a matrix and B\, B2, . ■ . , B a are 

b x b matrices. Let M be the ab x ab block matrix: 



M = 



/c 11 B 1 c 12 B 2 



\c al B 1 c a2 B 2 



ClaB a \ 



Then det(M) = ± det{A) b det(Bi) det(S 2 ) ■ ■ ■ det(B ). 

The matrix M is similar to a tensor product matrix, but the Bi may be different 
matrices, so this construction is somewhat more general. 



Proof. If det(A) = or det(-Bi) = for some i, then det(M) = as well. So, we 
assume that det (A) ^ and det(-Bj) ^ for all i. In order to find the determinant 
of M, we may transform M into a block upper triangular matrix using blockwisc 
row operations, obtaining a matrix M' in the following form: 



/c'n-Bi 




M' 



\ 



V 







C ' 22 B2 











CaaB, 



in which the c' u 's are the same as the entries obtained by the corresponding row 
operations applied to the matrix A. Now we have a block-upper triangular matrix. 
This implies that the determinant of the matrix is the product of the determinants 
of the diagonal block entries. Thus, 



h det(Bi) • (c' 22 ) b det(B 2 
) fc dct(Bi)det(B 2 ) 



det(M') = (c' n 

= ( c ll ^aa 

= det(i4') 6 det(Bi)det(B 2 )- 



■(0 6 det(B ) 
det(S„) 

det(B„) 
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We know that det(A') = ±det(A) (some row interchanges might have been neces- 
sary in the reduction to upper triangular form). So, we can substitute to yield: 

det(M) = ± dct(A) b det(-Bi) dct(B 2 ) • • • det(B a ), 

which is what we wanted to show. □ 

When S is a (k + 1) x {£ + 1) configuration, consisting of I + 1 distinct points 
on each of k + 1 distinct lines x = a u , u = 1, . . . ,k + 1, then V(P^^\ S) has the 
form given in the proposition, where A is the ordinary (k + 1) x (k + 1) univariate 
Vandermonde matrix of the x — a u , and B u is the (£ + 1) x + 1) univariate 
Vandermonde matrix of the y-coordinates of the points on x = a u . Hence we 
obtain the following consequence. 

Corollary 2.3. Let S be a (k + 1) x (£ + 1) configuration in (F*) 2 . Then 

detV(P k j;S) t^O. 

Proof. This follows from the factorization of the univariate Vandermonde determi- 
nants as products of differences of the x- and y-values at the points in S. □ 

We are now ready to prove the first major result of this section. 

Theorem 2.4. Let k,£ < q — 1 so that Pkj C Dq-i C M 2 . Then the minimum 
distance of the two-dimensional toric code Cp k l is 

d(C Pk J = (q- l) 2 - (k + £)(q -l) + k£= ((q 1) - k)((q 1) - £). 

Proof. We write d = d(Cp k e ). In order to show equality, we will show that both 
d< \q-l) 2 -(k + £)(q-l) + k£ &nd d> (q— l) 2 - (k + £)(q—l) + k£. We start with 
the former. In L — Span{l, x, x 2 , . . . , x , y, xy, x 2 y, . . . , x k y, . . . , x k y e } consider a 
polynomial p(x, y) = qi(x)q2(y) where qi(x) and qi{y) factor completely: 

p(x,y) = (x- ai)(x -a 2 )---{x- a k )(y - bi)(y - b 2 ) ■ ■ ■ {y - be), 

for some distinct a u and b v in ¥* This means that the codeword ev(p) has zeros 
in the positions corresponding to the points along k distinct vertical lines and £ 
distinct horizontal lines in (F*) 2 . We see that ev(p) has (k + £){q — 1) — k£ zeroes. 
Accordingly, we have d < (q - l) 2 - (k + £)(q - 1) + k£. 

Now we show d> (q— l) 2 — (k + £)(q — 1) + k£. By Corollary 2.3 and Proposition 
2.1, it suffices to show that every set T of size (k+£)(q—l) — k£+l in (F*) 2 contains 
a (k+l) x {£+!) configuration S. Since q~l > k, we see that (k+£)(q-l)-k£+l > 
£(q — 1). So, by the pigeonhole principle, some vertical line x = a\ contains £ + 1 
distinct points from T. There are at most (q—l) — (£+l) = q — £ — 2 other points 
on that line. Therefore, there are at least (k + £ — l)(q — 1) — k£ + 1 other points 
in T not on x = a\. To repeat this argument to find the rest of the configuration 
S, we must show that for all 1 < j < k, (k + £ - j)(q - 1) - k£ + 1 > £(q - 1 - j). 
To do this, we subtract £(q — 1 — j) and perform some arithmetic: 

(k + £- j)(q - 1) - k£ + 1 - £(q - 1 - j) = (k — j)((q -l)-£) + l>0. 

So, after we have found the first j sets of £ + 1 points on vertical lines in T, there 
are still enough additional points left in T so that we can find the next set of 
£ + 1 points. After k+l steps, we have a complete (k + l) x (£+ 1) configuration 
S c T. □ 
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The method used in the proof of Theorem 2.4 extends without difficulty to toric 
codes constructed from k\ x fc 2 x • • • x k m rectangular polytopes Pk lt ...,k m C W n 
for all m > 2. The result is as follows. 

Theorem 2.5. Let ki,...,k m be small enough so that Pk x ,....k m C Dg-i C M m . 
TTien £/ie minimum distance of the m- dimensional toric code Cp ki km is 

m 

d(C Pkl _ k J = l[((q-l)-k i ). 

i=i 

We next turn to the toric codes Cj^( m ) for Pe(m) an m-dimensional simplex of 
the form 

Pe{m) = conv{0, £ei, . . . , £e m }, 

where the ej are the standard basis vectors in R m . The monomials corresponding to 
the ( m ^~ £ ) integer lattice points in Pt{m) are all of the monomials in m variables of 
total degree < £. Naturally enough, the corresponding Vandermonde matrices arise 
in the study of multivariate polynomial interpolation using polynomials of bounded 
total degree. The next recursive definition gives the special configurations of points 
where we will be able to compute Vandermonde determinants for the Pe(m). 

Definition 2.6. Ifm = 1, anlth order simplicial configuration is any collection 
of distinct points in F*. For m > 2, we will say that a collection S of ( m ^ 1 ) 

points in (¥*) m is an m- dimensional Ith order simplicial configuration if the 

following conditions hold: 

(1) For some i, 1 < i < m, there are hyperplanes Xi — a\,X{ — 02, ... ,Xi = 
ae + i such that for each 1 < j < £ + 1, S contains exactly C"" 1 ^/ -1 ) points 
with Xi = dj. (Note that 

(m + A ^ (m - 1 + j - 1\ 

by a standard binomial coefficient identity.) 

(2) For each j, 1 < j < £+1, the points form an (m—X)- dimensional 
simplicial configuration of order j — 1 . 

We call these special configurations of points simplicial configurations because 
they mimic, to an extent, the arrangement of the integer lattice points in the cor- 
responding simplex Pi(m). For instance, Figure 1 shows a 2-dimensional simplicial 
configuration of order 2 in (Fjf ) 2 . ft consists of six points. (We write a for a prim- 
itive element in the field Fs.) Note that there is 1 = ( ~t ) point in S on the line 
x\ = a\ = a 4 , 2 = points on the line x\ — 02 = a 3 , and 3 = ( 1 ^ 2 ) on the line 

xi = 03 = a. 



(2) 
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Figure 1. A 2-dimensional simplicial configuration S of order 2. 



In order to state our next result, a sort of recurrence relation for the Van- 
dermonde determinants dct V(Pi (m); S) where S is a an m-dimensional simplicial 
configuration, we introduce some notation. 

Let S be an m-dimensional £th order simplicial configuration consisting of ( m ^~ ) 
points, in hyperplanes x rn = ai, . . . ,x rn = a^+i- Write S = S' U S" where S' is 
the union of the points in Xj = ai, ■ ■ ■ , at, and S" is the set of points in Xi — a(,+i- 
Also, let 7r : F™ — * F™ -1 be the projection on the first m — 1 coordinates. By the 
definition, it follows that both S' and tt(S") are themselves simplicial configurations, 
with S' of dimension m and order £—1, and tt(S") of dimension m — 1 and order £. 

Theorem 2.7. Let Pi(m) be as above and let S be an £th order simplicial config- 
uration of ( m ^ e ) points as in the paragraph above. Then writing p = (p%, . .. ,p m ) 
for points p £ (F*) m , 



det V(Pt(m); S) = ± JJ ( Pm - a i+1 ) det V(P t -i(m); S') det V(P e (m - l);7r(S")). 
pes 1 

Before we give the proof of this theorem, we will do two things. First, we give 
an example to illustrate what the theorem is saying. The idea for this computation 
comes from [J, where corresponding sets of points in R" 1 are identified as poised 
sets for interpolation by polynomials of degree bounded bounded by £. Consider 
all polynomials of degree < 2 in three variables and the Vandcrmonde matrix 
V(P2{3); S). For notational simplicity, write points in a 3-dimensional simplicial 
configuration S C (F*) 3 of order 2 as {xi,yi, z 4 ), for i = 1, . . . , 10 = ( 3 2 2 )- Here S' 
consists of the first four points in S, and 5"' consists of the other six points. Under 
the hypothesis that S is a simplicial configuration, we have z$ = z& = ■ ■ ■ = Z\q = c 
for some c = a^. Noting this, but ignoring other equalities between the coordinates, 
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we see V{P 2 (3);S) = 
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To evaluate the determinant of this matrix, we perform row operations to introduce 
zeroes. First subtract c times row 4 from row 10, then c times row 3 from row 9, 
c times row 2 from row 8, and finally c times row 1 from row 4. After rearranging 
rows, we find a matrix with a block of zeroes: 



(3) 



Zl — c 
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Up to a sign, the determinant of V(P2(3); S) is therefore equal to the product of 

/ zi - c z 2 — c z 3 — c z 4 — c \ 
xi(zi-c) x 2 {z 2 -c) x 2 {z 3 -c) X 4 {z 4 -c) 
yi(zi-c) 2/2(22-0) y 3 (z 3 -c) y 4 (z 4 -c) 

\zi(zi-c) z 2 (z 2 -c) z 3 (z 3 -c) Z 4 {Z4-C)J 

which equals 

(zi - c)(z 2 - c)(z 3 - c)(z 4 - c) dot V(Pi(3), S 1 ), 

and the determinant of the lower right 6x6 block in (3), which is V(P2{2), n(S")). 
Hence 

4 

detV(P 2 {3);S) =±Y[{ Zj -c) detV(P 1 {3),S')detV{P 2 (2),n{S")), 

as in the statement of the theorem. 

Second, we note the following immediate consequence of the theorem. 

Corollary 2.8. Let Pg{m) be as above and let S be an iih order simplicial config- 
uration of ( m + l ) points. Then det V{P t {m); S) ^ 0. 
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Proof. This is seen easily by a double induction on m and £. In the base cases 
m = 1, £ arbitrary, we have an ordinary univariate Vandcrmondc determinant, 
which is nonzero by the definition of a simplicial configuration. For the induction 
step, the recurrence given in the theorem then establishes the corollary since all the 
factors are nonzero under the hypothesis that S is a simplicial configuration. □ 

We will now give a proof of Theorem 2.7, following the methods from 0j used 
in the example above. 

Proof. Since S is simplicial, the ( m ^f 1 _1 ) points in S" C S all have the same 
x m -coordinate, say x m = a^+i = c. For each pair of monomials x e and x e x m 
corresponding to points in Pi(m), we perform a row operation on the Vandermonde 
matrix, subtracting c times the row of V(Pi(m), S) for x e from the row for x e x rn in 
decreasing order by the degree in x m . After rearranging rows to put all the zeroes 
created by these operations in the upper right block, the lower right block in the 
columns corresponding to S" is the matrix V(Pe(m — l),ir(S")). In the upper left 
block (in the columns corresponding to 5"), all the entries in a column are divisible 
by one of the (p m — c) for p E S'. Factoring out those factors from each column, the 
matrix that is left is V(Pe-i(m), S'), and the statement of the theorem follows. □ 

We note that if we apply the recurrence relation repeatedly, another corollary of 
Theorem 2.7 is a closed formula for det V(Pz(m); S) for S a simplicial configuration 
in terms of univariate Vandermonde determinants. We will not need this formula, 
so we leave the derivation of its exact form as an exercise. 

We will now use Corollary 2.8 to establish the minimum distances of the toric 
codes C Pt{m ). 

Theorem 2.9. Let I < q — 1, and let Pi(m) be the simplex in R m defined above. 
Then the minimum distance of the toric code Cp t r m ) is given by 



inequalities d < (q - l) m - £(q - l)™" 1 and d > (q - l) m - £(q - l)" 1 " 1 hold. The 
first of these follows as in Theorem 2.4, since the completely reducible polynomials 



for aj distinct in F* arc contained in L — Spanjx 6 : e € Pf (m)}. Such a polynomial 
has zeroes at all the F g -rational points with nonzero coordinates on the union of 
the hyperplanes x m = aj. There are £(q — l) m_1 such points, so d < (q — l) m — 
£(q — l)™ -1 as claimed. 

To establish the reverse inequality, by Corollary 2.8 and Proposition 2.1, it suf- 
fices to show that every set T of £(q — l) m_1 + 1 points contains an m-dimensional 
simplicial configuration S of order £. By the pigeonhole principle, T contains some 
set of £(q — l)™- 2 + 1 points with the same x m -coordinate x m — ae+i. By an easy 
induction on m, it follows that 



(the first term on the right of (2)) for all m > 2 (and q — 1 > £, of course) There are 
at least (£— + 1 points in T that do not lie on x m = a^+i. Hence we can 
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apply the same pigeonhole principle argument repeatedly (or argue by induction 
on £) to see that T contains 

g( m ";- + r i )-( m « + o 

points making up an m-dimensional £th order simplicial configuration. □ 

We also have the following consequence for more general simplices. Let £ 1 > 1 
for all i and define 

Pe 1 ,...,t m = conv{0,£iei, . . . ,£ m e m }, 
where again the ej are the standard basis vectors in R m . 

Corollary 2.10. If Pi u ...,i m C CM™, and £ = mas, l u then 

d(C Pei _ e J = (q-ir-£(q-ir- 1 . 

Proof. By definition, Pg 1 ....j m Q Pt(m). Hence Cp, is a subcode of Cp f ( m ) and 

d(C Pei _ e J > d(C Pdm) ) = {q-l) m -£(q-l) m -\ But if £ = £ io , then C Pii _ 4m 
also contains codewords of weight exactly (q — l) m — £(q — l) m_1 obtained from 
evaluation of completely reducible polynomials 

p(xi) = (x io - on) • • • (x la - a t ) 

for some i and distinct a,j € F*. This establishes the corollary. □ 

The same sort of reasoning applies to any toric code from a polytope P C Pe(m) 
that contains one complete edge conv{0,tei} or conv{£ej, £ej} of the simplex P C 
Pt(m) and gives the same minimum distance for Cp. 

3. Classification of toric codes 

In this section we will begin by stating and proving a theorem guaranteeing 
that two toric codes have the same parameters. We begin by introducing some 
terminology 

Definition 3.1. Let C\ and C2 be two codes of block length n and dimension k 
over F q . Let G\ be a generator matrix for C\. Then C\ and C'2 are said to be 
monomially equivalent if there is an invertible n x n diagonal matrix A and an 
n x n permutation matrix II such that 

G 2 = GiAU 

is a generator matrix for C2 ■ 

It is easy to see that monomial equivalence is actually an equivalence relation 
on codes since a product LTA equals ATI for another invertible diagonal matrix A'. 
It is also a direct consequence of the definition that monomially equivalent codes 
C\ and C 2 have the same dimension and the same minimum distance (indeed, the 
same full weight enumerator). 

Next we turn to a natural notion of equivalence for polytopes. Recall that an 
affine transformation of R m is a mapping of the form T(x) = Mx + A, where A is 
a fixed vector and M is an m x m matrix. The affine mappings T where M, A have 
integer entries and M e GL(m, Z) (so det(M) = ±1) are precisely the bijective 
affine mappings from the integer lattice Z m to itself. 
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Definition 3.2. We will say that two integral convex polytopes P\ and P2 in R m 
are are lattice equivalent if there exists an invertible integer affine transformation 
T as above such that T(P\) = P2. 

This brings us to our next theorem which relates the two concepts we have just 
defined. 

Theorem 3.3. If two polytopes P\ and P2 are lattice equivalent, then the toric 
codes Cp 1 and Cp 2 are monomially equivalent. 

Proof. Suppose we have two lattice equivalent polytopes P\ and P2. Both Pi and 
P2 contain integer lattice points corresponding to monomials of the form x e where 
e £ Z m . By our hypothesis on Pi and P2, there exists an invertible integer trans- 
formation 

T(x) = M(x) + A 

such that T(P 1 ) = P 2 and M is an element of GL(m,Z) so det(M) = ±1. 
Hence #(Pi) = #(P 2 ). Let P 1 n Z m = {e(i) : i = 1, . . . , #(Pi)}. So, C Pl is 
spanned by ev(x e ^) for 1 < i < n, and similarly Cp 2 is spanned by ev(x T ^ 1 ^). 
Write a for a primitive element in ¥*. Let e(i) € Pi n Z m and define a? — 
(a? 1 , . . . ,a^ m ) e (F*) m . The component of ev(x e ^) e C Pl corresponding to a? is 
a ( e M>/} ; where (e(i),f) is the usual dot product. The corresponding entry in the 
codeword ev(x T ^ 1 ^) in Cp 2 is is a( T ( e W)>/> . This can be rewritten as 

a (Me(i) + \,f) = a (Me(i),f) . Q (A,/) 

The second term of the product is not dependent on e(i). These nonzero scalars are 
the diagonal entries in the matrix A as in the definition of monomially equivalent 
codes. By a standard property of dot products, 

a (Me(i)J) = a (e(i),M*/>_ 

The transposed matrix M* also defines a bijective mapping from Z m to Z m since 
det(M') = det(M) = ±1. Now we must show that M* induces a permutation of 
(¥*) m . Suppose M l f ee M l g (mod q - 1). Since det(M*) =±1^0, we know 
that M* is invertible and (M*) -1 is also an integer matrix. So, we can multiply by 
(M*) _1 on the left. Hence, / = g (mod q — 1) and M t defines a permutation of 
the points a? , as desired. Note that M* permutes all of the codewords in the same 
way. This gives the permutation matrix n. Hence Cp 1 is monomially equivalent to 
C P2 . □ 

In the remainder of this section, we will show how this result leads to a complete 
classification for toric codes with m = 2 and k < 5. The classification could also be 
continued, of course, using a census of lattice equivalence classes of lattice polytopes 
with given #(P). 

Proposition 3.4. Every toric surface code Cp with k = 2 is monomially equivalent 
to the toric code Cp 2 for P2 = conv{(0, 0), (1,0)}. 

Proof. Let e(l), e(2) G 1? be the integer lattice points in P. We can use a transla- 
tion to map e(l) to (0, 0). Then let e(2) = (a, b) e J?. By convexity we have that 
gcd(a, b) = 1, since otherwise there would be additional integer lattice points on the 
line from e(l) to e(2) and k = #(P) would be greater than 2. Since gcd(a, b) = 1, 
there exist integers r, s such that ra + sb= 1 , and this implies that there exists an 
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invertible integer matrix M = 



such that M 



. Hence there 



-b a J \bj V°y 

is an affine equivalence between P and P 2 . By Theorem 3.1, this completes the 
proof. □ 

Next, we wish to find a "nice" lattice polygon in each possible lattice equivalence 
class with = 3, 4, 5. One way is to add additional points to P 2 - Using Pick's 

Theorem: A(P) = + \d(P) — 1, (where d(P) is the number of lattice points 

in the boundary of P) then eliminating lattice equivalent polygons, we obtain the 
following. 

Theorem 3.5. Every toric surface code with 3 < k < 5 is monomially equivalent 
to one constructed from one of the following polygons. 



V - 



>(i) 




XX 2 1 X X 

Figure 2. Polygons yielding toric codes with k = 3. 



,(2) 





Figure 3. Polygons yielding toric codes with k = 4. 
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Figure 4. Polygons yielding toric codes with k = 5. 

The final step in our classification is to show that no two of the toric surface 
codes constructed from these polygons can be monomially equivalent, hence they 
lie in distinct monomial equivalence classes. We do this by applying results from §1 
to show that the minimum distances (or in some cases, other parts of the complete 
weight enumerators) are distinct. 

Theorem 3.6. Let q > 5. No two of the toric codes Cp(¥ q ) constructed from the 
polygons in Theorem 3.5 are monomially equivalent. 

Proof. If the dimensions are different, the toric codes are certainly not monomially 
equivalent. Hence we only need to consider each k separately. 

For the two codes with k = 3, Theorem 2.4 (which applies when I — also) 
shows d(C p (i)) = (q — l) 2 — 2(q — 1). On the other hand, Theorem 2.9 gives 

d(C p <2)) = (q — l) 2 — (q — 1). Hence these two codes are not equivalent. 

For the codes with k = 4, Theorem 2.4 shows d(C pW ) = (q - l) 2 - 3(q - 1). 
Corollary 2.10 shows d(C p <2)) = (q — l) 2 — 2(q — 1). Theorem 2.4 applies to C p (?,) 
also, and shows d(C < 3 )) = {{q — 1) — l) 2 = (q — l) 2 — {2q — 3). Finally, we must 
analyze d(C <4>). Write C (4)(F 9 ) = C(¥ q ). In this case, some more advanced 

4 ' 4 

tools are needed. If we translate this polygon by (1, 1) to place it in then 
we are evaluating polynomials in Span{l, xy, x 2 y, xy 2 } to get the codewords of 
the corresponding (monomially equivalent) code. Any linear combination of these 
monomials in which the coefficient of x 2 y or xy 2 is nonzero defines an absolutely 
irreducible curve of degree 3, whose closure in P 2 has arithmetic genus 1 by Theorem 
4.2 of Hence by the general version of the Hasse-Weil bound from p^, there can 
be at most 1 + q + 2^/q F g -rational points on the corresponding affine curve. This 
means that the minimum distance of C is at least (g — l) 2 — (1 + q + 2^/q). On 
the other hand, the other k = 4 examples have minimum distance no larger than 
d(C „(3)) = (q — l) 2 — (2q — 3). It is easy to see from the quadratic formula that 
(1 + q + 2^/q) < 2q — 3 for all q > 11. Hence d(C) is strictly larger than any of the 
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others for q > 1 1 . For the remaining small values of q we check directly that d(C) 
is different from the others using the Magma code described in [H] . The results are: 



9 
11 



d(C(F,)) 



27 
40 
52 
85 



(q-l) 2 -(2q-3) 



25 
3G 
49 
81 



Table 1. d(C pW (F,)) 



These are also different from any of the other k — 4 codes over those fields. 

For the k = 5 codes, Theorem 2.4 shows d(C p m ) = (q - l) 2 — 4(g— 1). Corollary 



2.10 shows d(C p (2) 



(q 



3(9 



1). C p (3) 



is a subcode of the code C p <2) 



from Theorem 2.9, and contains codewords of the same minimum weight as the 
supercode. So d{C p <s) ) = (q— l) 2 — 2(q — 1). The other three k = 5 codes also have 

d(C p (i)) — (q — l) 2 — 2(q — 1), which can be seen, for example, using MinkoWSki- 
^5 

decomposable subpolytopes (the sets of three collinear points) as in ^U]. In C p (4>, 
for example, we have codewords ev(b(xy — a\){xy — 02)), where ai,a,2,b € F* and 
a\ ^ 0,2, which have 2(q — 1) zeroes in (F*) 2 . 

To show that the four codes with d = (q — l) 2 — 2(q — 1) are not equivalent, 
we need to look at finer invariants. For instance, C p (b) can be distinguished from 

(5) 

the other three by the number of words of minimum weight. In P 5 , there are two 
different sets of three collinear lattice points while in the others, there is only one. 
This means that there will be more words of the minimum weight in C p (a) than 



in C_(i) for i 



3,4,6. C„(5) 



has at least 2( q 2 ls j(q — 1) such words because there 
are two distinct families of reducible polynomials: b(x — a\)(x — a^) with b, a; € ¥* 
and ai 0,2 and b(y — ai)(y^ 1 — a 2 ) b, cii € F* and ai 7^ a 2 _1 . On the other 
hand, C p (%) for i = 3, 4, 6 have (at least) ( l? 2 1 )(o , ~ -0 sucn words. There are more 

for some small q, but never as many as 2( q ~ 1 )(q — 1). See the weight enumerators 
for C p {6) over Fn and F16 in Table 2 below. For sufficiently large q, we claim in 

fact that there are exactly ( 9 ^ 1 ){<1 — 1) such words. This follows from the general 
Hasse-Weil bound from pQ. For instance, for C p ( 6 ) , if q is sufficiently large, then 

we claim all words in C„(6) of weight (q — l) 2 — 2(q — 1) come from evaluations 

ev(b(y — ai)(y — ^2))- Any other such word could come only from evaluating a 
linear combination of {1, x, y, y 2 , x -1 ?/ -1 } in which y 2 , x, x~ l y~~ l all appear with 
nonzero coefficients (since otherwise we are in a case previously covered). Any 
such curve is absolutely irreducible, of arithmetic genus 2 (because of the 2 interior 
lattice points in this case, see Theorem 4.2 of 0). A simple argument shows that 
1 + q + 4^/q < 2q — 2 for all q > 23. For smaller values of q, we verify directly that 
the weight enumerators of C p <5) (F q ) do not match the weight enumerators of the 

other codes using the Magma code from See Table 2. 

To distinguish C p <3) and C p <4) , we use the codewords of weight (q— l) 2 — (2q — 3) 

5 5 

(one more than the minimum weight). Both of these codes contain such words 
coming from evaluation of the polynomials corresponding to the lxl squares 
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contained in these polygons (copies of P4 ). Any such square yields (q — l) 3 words 
of this weight since the polynomials in question have the form c(x — a)(y — b) and 
a, 6, c S V* are arbitrary. However C p <4) has precisely (q — l) 3 words of weight 

(q — l) 2 — (2q — 3), while C p (3> has more of them, 3(o — l) 3 to be specific. This can 
be seen by considering the reducible polynomials d(x — a)(y — bx — c) that evaluate 
to give codewords in C„( 3 ) . We get codewords of weight (q — l) 2 — (2q — 3) if b = 0, 
or if c = 0, or if 6, c 7^ and a = —c/b. 

Finally, to distinguish C p (e) from the other three codes with d = (q— l) 2 — 2(q— 1), 
we must argue as in the last case of the k — 4 codes. If q is sufficiently large, then we 
claim C p <6) contains no words at all of weight (q — l) 2 — (2q — 3) . By Corollary 2.10 
and the previous cases, we see that any such word could come only from a linear 
combination of {1, x, y, y 2 , x~ 1 y~ 1 } in which y 2 , x, x~ 1 y~ 1 all appear with nonzero 
coefficients. As before, any such curve is absolutely irreducible, of arithmetic genus 
2. A simple argument shows that 1 + q + 4^/q < 2q — 3 for all q > 23. Hence the 
Hasse-Weil bound from t h shows that there are no words of this weight for large 
q. For smaller values of q, we again verify directly that the weight enumerators of 
C p (6) (¥ q ) do not match the weight enumerators of the other codes. See Table 2 
below. 

The following table gives the first three nonzero terms in the weight enumerators: 

(«-i) 2 
W c (x)= ]T A z x\ 

i=0 

where Ai = \{w G C : wt(w) = i}, for the k — 5 toric codes with d — (q — l) 2 — 
2 (q — 1). These were all computed using the Magma code from [Hj. 

Over F 7 : 

1 + 90x 24 + 648x 25 + • • • 
1 + 90x 24 + 216x 25 + ■■■ 
1 + 18Cte 24 + 324x 26 + • • • 
1 + 90a; 24 + 432a; 26 + ■ • • 



1 + 147a; 35 + 1029a; 36 H 

1 + 147.x 35 + 343a; 36 H 

1 + 294a; 35 + 343a; 37 H 

1 + 147a; 35 + 1029a; 37 H 



1 + 224a; 48 + 1536a; 49 + • • ■ 
1 + 224.x 48 + 5 1 2x 49 + • • • 
1 + 448x 48 + 5 1 2x 51 + • • • 
1 + 224x 48 + 5 1 2x 50 + • • • 



(3) 
U) 
(5) 
(6) 



Over 



,(3) 
5 

,(4) 
5 

,(5) 
5 

,(6) 



Over Fq 



p(3) 

^5 

p(4) 
p(5) 



P. 



(6) 
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Over Fn 



Over F13: 



Over Fi 6 : 



Over F17: 



Over Fi S 



Over F23: 



P 5 (3) 1 + 450a; 80 + 30 00.x 81 + • 

P 5 (4) 1 + 450a; 80 + 1000a; 81 + • 

P 5 (5) 1 + 900a; 80 + 1500a; 84 + • 

P 5 (6) 1 + 650a; 80 + 1000a; 82 + • 

P 5 (3) 1 + 792a; 120 + 5184a; 121 + ■ 

P 5 (4) 1 + 792a; 120 + 1728a; 121 + ■ 

P. (5) 1 + 1584a; 120 + 77 76a; 126 + 



P, (6) 1 + 792a; 120 + 1728a; 125 



P 5 (3) 1 + 1575a; 195 + 10125a; 196 + • • 

P 5 (4) 1 + 1575a; 195 + 33 75a; 196 + • • 

P 5 (5) 1 + 3150a; 195 + 13 5 00a; 203 + • • 

P 5 (6) 1 + 2250a; 195 + 13 5 00a; 203 + • • 

P 5 (3) 1 + 1920a; 224 + 12288a; 225 + • • 

P 5 (4) 1 + 1920a; 224 + 40 96a; 225 + • • 

P, (5) 1 + 3840a; 224 + 5120a; 232 + • • 



.230 



P, (6) 1 + 1920a; 224 + 40 96a; 



P 5 (3) 1 + 2754a; 288 + 17496a; 289 + • 
P, (4) 1 + 2754a; 288 + 58 32x 289 + • • 



P 5 (5) 1 + 5508a; 288 + 32076a; 298 
P, (6) 1 + 2754a; 288 + 58 32a; 294 - 



P 5 (3) 1 + 5082a; 440 + 3 1 944a; 441 + • ■ 

P 5 (4) 1 + 5082a; 440 + 106 48a; 441 + • ■ 

P 5 (5) 1 + 10164a; 440 + 154396a; 454 + 

3(6) 1 J- 5082a; 440 + 21296a; 450 ->- ■ 



5 



Table 2. 

Hence the enumerators never coincide for these four codes, even in exceptional 
cases for small q. □ 

References 

[1] Y.Aubry, M.Perret, A Weil theorem for singular curves, in Arithmetic, Geometry, and Coding 
Theory, R.Pellikaan, M.Perret, S.G.Vladut, eds. de Gruyter, Berlin, 1996, 1-7. 

[2] P. Beelen, R. Pellikaan, The Newton polygon of plane curves with many rational points, Des., 
Codes and Cryptography 21 (2000), 41-67. 

[3] A.E. Brouwer, Bounds on linear codes, Handbook of Coding Theory, Elsevier (1998), 295-461, 
and online updates at http://www.win.tue.nl/-aeb/voorlincod.html. 



ON m-DIMENSIONAL TORIC CODES 



17 



[4] C.K. Chui, M.J. Lai, Vandermonde determinant and Lagrange interpolation in R s , in Non- 
linear and convex analysis (Santa Barbara, Calif., 1985), 23-35, Lecture Notes in Pure and 
Appl. Math., 107, Dekker, New York, 1987. 

[5] V. Diaz, C. Guevara, M. Vath, Codes from n- Dimensional Polyhedra and n-Dimensional 
Cyclic Codes, Proceedings of SIMU summer institute, 2001. 

[6] M. Gasca, T. Sauer, Polynomial interpolation in several variables. Multivariate polynomial 
interpolation., Adv. Comput. Math. 12 (2000), 377-410. 

[7] J. Hansen, Toric surfaces and error- correcting codes, in Coding theory, cryptography and 
related areas (Guanajuato, 1998), 132-142, Springer, Berlin, 2000. 

[8] J. Hansen, Toric varieties Hirzebruch surfaces and error- correcting codes, Appl. Algebra 
Engrg. Comm. Comput. 13 (2002), 289-300. 

[9] D. Joyner, Toric codes over finite fields, Appl. Algebra Engrg. Comm. Comput. 15 (2004), 
63-79. 

[10] J. Little, H. Schenck, Toric surface codes and Minkowski sums, preprint, 2005. 
[11] G. Ziegler, Lectures on Polytopes, Springer Verlag, Berlin, 1995. 

Department of Mathematics and Computer Science, College of the Holy Cross, 
Worcester, MA 01610 

E-ma il address: |lit tle@mathcs .holycross . ed u| 

URL: http: //mathcs .holycross . edu/~little/homepage . html 

Department of Mathematics and Computer Science, College of the Holy Cross, 
Worcester, MA 01610 



